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Abstract 

In this note, we find sufficient conditions for an operator with kernel of the form 
A(x)B(y) — A(x)B(y)/(x — y) (which we call a Tracy- Widom type operator) to be the square of a Hankel 
operator. We consider two contexts: infinite matrices on I 2 , and integral operators on the Hardy space 
H 2 (T). The results can be applied to the discrete Bessel kernel, which is significant in random matrix 
theory. 
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1 Introduction 

In random matrix theory it is natural (see, e.g. [I]) to consider integrable operators T, where the kernel of 
T is 

A j (z)B j (w) 



3=1 

and Y^j=i Aj(z)Bj(z) = 0. Here we are concerned with a special class of such operators, namely those with 
kernel of the form 

x-y 

which we shall refer to as Tracy-Widom operators. The variables x and y may be non-negative integers, as 
in the discrete kernels considered in section 2, continuous real parameters, as in e.g. [2], or may live on the 
circle, as in section 3. We look for conditions under which these operators can be expressed as T 2 or T*r, 
where T is a Hankel operator. In particular we recover a result of Borodin et al [3], showing that the discrete 
Bessel kernel can be written as 



x — y 

y fc=0 



(3) 



We can then read off information about K from knowledge of the Hankel operator T. For example, a trace 
formula follows immediately, and the spectrum of K can be calculated from the spectrum of T (which in 
many cases is easier to calculate). Megretski, Peller and Treil [I] have characterised the self-adjoint bounded 
linear operators that are unitarily equivalent to Hankel operators: we apply their results to gain spectral 
information about the operators K. 



2 Discrete integrable operators 



Define No = N U {0}. We consider infinite matrices with kernel K(x,y), where K(x,y) is defined by @. 
Recall that a Hankel matrix = [<f>(m + n)] m . n >o with (</>(fc)) € £ 2 has square 



<^>(m + k)4>(n + k) 



.k=0 



(4) 



Nehari's theorem (see, e.g. [5j p. 3]) states that is a bounded operator on € 2 (No) if and only if (4>(n)) are 
the positive Fourier coefficients of some function in L°°(T). We write the kernel K(x,y) in matricial form, 

1 



K(x,y) 
a(x) 



x-y 



(Fa(z),a(2/)>> (x^y) 



A(x) 


, F = 


-1 




. B(x) . 




1 



(5) 



and look for sufficient conditions under which we can construct a function <j> : No — > C with (4>(j)) € Z 2 , such 
that 

oo 

K(x,y) = ^</>(x + k)<l>{y + k), (x^y). (6) 



fe=0 



Definition 2.1 Lef 5 be the shift operator on ^ 2 (No), so that Sf(x) = f(x—l) (where we define /(— 1) = 0), 
and let R be the adjoint shift operator Rf(x) — f(x + 1). The forward difference operator A is defined by 
Af(x) = f(x + 1) — f(x). Notice that A = R x — I. Where there are several variables, we write R x , A y and 
so on. 

As usual, A T is the transpose of a matrix A, while B* denotes the adjoint of an operator B. 
Lemma 2.2 (Lyapunov equation) Suppose that R and B are bounded linear operators on i 2 such that 

oo 

(R J BB*(R*) j £,,^) < oo for all £ e £ 2 , 

so that the series 



3=0 



K = ^R>BB*(R*y 

j=o 

is convergent in the weak operator topology. Then 

K — RKR* = —BB*. (7) 

Proof. Clear from calculation of the left hand side of (0 • ■ 

In the following Lemma, we state explicitly the specialisation of the above result to discrete kernels. 

Lemma 2.3 Let &(x,y) be any function $ : N 2 , — > C, and suppose <fr : No — * C is such that (<fi(j)) £ Q 2 ■ 
Then 



<&(x, y) = <j)(x + k)<f)(y + k) for all x, y e N 



(8) 



fe=0 



if and only if 

(A x S y + A y )$(x, y) = -<j>{x)(j){y) for all x, y e N (9) 

and 

&(x,y)^0 as x or y — ► oo. (10) 
Proof. Suppose (O holds. Then we have 

oo oo 

(A x S y + Ay) 4>(x + k)4>(y + k) = (S x S y - I) ^ <p(x + k)(f>(y + k) 

k=0 k=0 

oo 

= (<1>( X + k + l )^y + k + l)-(t)(x + k)(f>(y + k)) 

k=0 

= -<f>(x)<t>(y), (11) 
so that © holds. By the Cauchy-Schwarz inequality, and since (</>(j)) £ £ 2 , we have 

oo 

$(x,y) = ^(j}{x + k)(t){y + k) 

k=0 

/ oo \ Va / oo \ V2 

< ^(x + fc) 2 J>(y + fc) 2 -0 as zory^oc, (12) 

\fc=0 / \fe=0 / 

which is condition (fT0|) . Conversely, suppose that we have ([9]) and (fTX)|) . and let 

oo 

G(z, y) = $(2;, y) - £ 0(x + fc)0(y + fc). (13) 

fc=0 

By ®, we have 

(A x S y + Ay)G(x,y) = for all x, y S N , 

so that G(x, y) = G(x + 1, y + 1) for all 1,1/ £ Nq. We then use the hypothesis (jTUJ) and the estimate in (fT2"]) 
to show that G(:r, y) — > as a; or y — > 00, and hence that G is identically zero for all non-negative integers 
a; and y, so that ([8]) holds. ■ 

Theorem 2.4 Let K(x,y) be as defined in with (a(x))!^ ~ ([A(x), B(x)] t ) x < Lq a sequence of 2 x 1 
real vectors such that 

]T||a(*)|| 2 <co. (14) 

a;>0 

Suppose that there exists a sequence of 2 x 2 reaZ matrices S x such that &(x + 1) = S x &(x) for all x £ No 
and that 

S T FS X - F 

(15) 

x-y 

is a constant matrix. Then C is symmetric. Suppose further that C has eigenvalues A € IR \ {0} and 0, and 
let [a,/?] T be a real unit eigenvector corresponding to A. Then 

00 

K(x,y) = -sgn (A) ^ <f>(x + k)4>(y + k) for x,yeN (x ^ y), (16) 

fc=0 

where 

<P{x) = \\\ 1/2 (aA{x)+f3B{x)) (17) 



and {(j>{x)) e i 2 



Proof. We set 



S„ FS X — F 
c = _y * 

x-y 



(18) 



where C is constant by hypothesis, so that we can exchange the roles of x and y, and find that C T = C. We 
have, for i^y, 

(A x S y + A y ) K(x,y) = (S x S y -I)^—(Fa(x),a(y)) 

x y 

= S x l — (F a (x),S ya (y)) - -J— (Fa(i), a(y)) 

x — y — 1 x — y 

= — (FS x a(x), S„a(y)) - — (Fa(x), a(y)> 
a; — y x — y 



= ^((^FS x -F)a(x),a(j/)) 
x — y 

= (Ca(x)Mv)) ■ 



(19) 



Since C is real and symmetric, and by hypothesis has eigenvalues A ^ and 0, there exists a real orthogonal 
matrix U of unit eigenvectors such that 



U 1 CU = 



A 




We have 



(20) 



(A x S v + A v )K(x,y) = (Ca(aO,a(y)> 

A 



U 





1 




A 



1 


sgn(A)</>(x)</>(2/), 



C/ T a(x),a(y)\ 
E/ T a(x),t/ T a(<,) 
£/ T a(x), 



1 




U T a(y) 



(21) 



where 



4>(x) 




|A| 1/2 











J7 T a(x). 



(22) 



Note that (4>{x)) S by the condition X^>o ll a ( x )l| 2 < °°> since [/ is a constant matrix. It is also clear that 
K(x,y) — > as a; or y — > oo, by the same condition on a(x). We now let [a,/3] T be a real unit eigenvector 
of C corresponding to A, and the result follows by Lemma l2~3l ■ 



Corollary 2.5 Let K(x,y) be as defined in {5p, with (a(x)) x < L = (L4(x), B(x)) T ) x < L a sequence of 2 x 1 
real vectors such that 



< oo. 



(23) 



x>0 



Suppose that a(x + 1) = {Lx + M)a(x) (for all x € NoJ, where L and M are real constant 2x2 matrices 
that satisfy 

det L = 
det M = 1 

M T FL is symmetric, and has eigenvalues A S IR \ {0} and 0. 
Lei [a,/3] T fee a real unit eigenvector of M T FL corresponding to A. Then 



K{x,y) = -sgn(X) ^2 <j)(x + k)(/)(y + k) for all x,y G N (x^y), 



fc=0 



w/iere = |A| 1/2 (aA(x) + (3B(x)), and (0(x)) € 



Proof. We have M T FM = F det M (indeed, this is true for any 2x2 matrix) and hence 
M T FM = F. Likewise L T FL = 0. Setting S x = Lx + M as in Theorem |2"H we now have 

S?FS X - F (Ly + M) T F(Lx + M) - F 



x-y 



x-y 

M T FLx - (M T FL) T y 



x-y 
M T FL(x - y) 



(since F T = -F) 



(since M T FL is symmetric by hypothesis) 



= M 1 FL. 



(24) 



(25) 



Hence C — (SyFS x — F)/(x — y) is a constant matrix. Thus, together with the summability criterion on 
the sequence (a(x)), the hypotheses of Theorem 12.41 are all satisfied, so we have the result. ■ 



Example 2.6 



Let J y {z) be the Bessel functions of the hrst kind of order i>, and write J x = J x (2\^9), where 9 is a 
positive real parameter. The discrete Bessel kernel 

J(x, y; 9) = Vd JxJy+1 - JyJx+1 



x-y 



(26) 



arises in the study of various discrete- variable random matrix models, as in [6] and [3]. Note that 
J x is an entire function of x, so that J(x, X] 9) is well-defined via L'Hopital's rule. In the notation of 
Corollary |2.5l we take 

" V9J X 

Jx+l 



a(x) = 



(27) 



The standard formula (see [U p. 379]) 

oo oc 

e t2tsine = J (2t) + 2 Jam (2*) cos 2m0 + 2i ^ J 2m _i(2t) sin(2m - 1)6* (28) 

m—l m— 1 

and Parseval's identity can be used to show that Jo(2t) 2 + 2 Ylm=i Jm(2t) 2 — 1 for all real t, and hence 
that the sequence (J x )%Lo is square summable. Thus the condition J2 x >o \\ a ( x )\\ < 00 is satisfied. 
The 3-term recurrence relation for the Bessel functions 



J x+2 (2z) - ?-tlj x+1 (2z) + J x (2z) = 



(29) 



becomes 

a(x + l) = 

and so we have a(x + 1) = [Lx + M)a(x), where 

L = 










Ve ' 




Vej x 


Jx+2 




-1 


x + 1 




Jx+l 



(30) 



and 



M 




° 

o Ve 

-i i 



It is clear that these matrices satisfy detL = and detM = 1, and we have 



M T FL = 





-1 



so we pick the unit eigenvector [a, 0\ T = [0, 1] T . Thus, the function </>(x) in Corollary [23] is J^+i (2\/#) , 
and we recover a result of Borodin et al in [3] 



J(x,y;9) = J2j x+k+1 (2V6)J y+k+1 (2Ve), x,y e No, 



(31) 



fc=0 



without their use of asymptotic formulae for the Bessel functions. 



The preceding results are identities of kernels for x ^ y. Evidently, the sum in the right-hand side of (|16[) 
makes sense for x = y, and hence gives one possible extension of the left-hand side to the case x = y. We 
use the extension to define an operator K with matrix given by K(x, y). 

Proposition 2.7 Suppose that the vector [A(x), B(x)) T satisfies the conditions of Theorem \2.4\ so that 
K(x,y) — X^fcLo fy<f*(y + k). Suppose also that ^^Lq n<fi(n) 2 < oo. Then the operator K represented 

by the matrix [K(x, 2/)]^ =0 is trace class and has trace: 



traced = ^(x + l)<f)(x) 2 



(32) 



x=0 



Proof. The summability condition on <f> ensures that T<p is Hilbert-Schmidt, which implies that K = T| is 
trace-class. We have 

oo oo oo 

trace K = ^ K(x, x) = ^T / ^T / <f>(x + kf (33) 



x=0 



x=0 k=0 



from which the result follows immediately. 



Definition 2.8 For a compact and self-adjoint operator W on a Hilbert space H , the spectral multiplicity 
function i/^y(A) : IR — > {0, 1, . . .} U {oo} is given by 



uw(X) = dim{x e H : Wx = Xx} (A € IR). 



(34) 



We now give the consequences of a result of Peller, Megretski and Treil in [J] in the case of discrete integrablc 
operators. 



Proposition 2.9 Suppose that T$ and K are as in Proposition \2.7\ Then and K are compact and 
self- adjoint, and 

(i) Vk(0) — or Vk(0) = oo; 

(ii) for A > 0, vk{X) < oo and vk{X) = J-T^v'A) + Vr+(— y/X); 
(Hi) if vk{X) is even, then z/r^VX) = vr^—VX); 

(iv) if vr{X) is odd, then vr^yX) — vy^ (— vA) = 1. 

Proof, (i) follows from Beurling's theorem (see [S], page 15), while (ii) is elementary. Peller, Megretski and 
Treil show in 4J that for any compact and self-adjoint Hankel operator T, the spectral multiplicity function 
satisfies ^(A) — ur(— A)| < 1. Using this, and (ii), statements (Hi) and (iv) follow immediately. ■ 

Remark 2.10 The Carleman operator V : L 2 (0,oo) — > L 2 (0, oo) is given by 

Tf(x)= [°° J— f(t)dt, (35) 
x + 1 



so r has kernel of Tracy- Widom type 



r 2 /H= r^^/(*)dt. (36) 

u — t 



Carleman showed that V is a positive self-adjoint Hankel operator with continuous spectrum [0, n] of multi- 
plicity two (see [5, p. 44%])> so ^ e Tracy-Widom type operator T 2 has spectrum [0, it 2 }, also of multiplicity 
two. This contrasts with (iii) and (iv) of Proposition[ 



3 Integrable operators on H 2 

Let H 2 be the usual Hardy space on the unit circle T, with orthonormal basis {1, z, z 2 , . . .}, and let R+ : L 2 — > 
H 2 and i?_ : L 2 — > L 2 H 2 be the Riesz orthogonal projection operators. We let denote multiplication 
by 4>, and define the Toeplitz operator on H 2 with symbol <fi to be T<j, = R+M$R+. Let J : L 2 — > L 2 be a 
flip operator, whose operation on a function / 6 i/ 2 is Jf(z) — zf(z). Note that J maps iJ 2 onto L 2 Q H 2 
(and vice versa) and that J 2 = I. The Hankel operator on H 2 with symbol € L°° is then 

T = Ji?_M . (37) 

We let the integral operator W on L 2 (T) have kernel 

.(.V^ g^^ , (38) 

where operates on a function / e L 2 (T) in the usual way: 

Wf(e ie ) = ^- [ W{e w ,e^)f(e^)^. (39) 



Lemma 3.1 Suppose that f,g 6 L°° have f = g. Then W defines a bounded and self-adjoint operator on 
L 2 . Further, R+WR+ : H 2 -> H 2 satisfies 

R+WR+ = r;r / -r;r g . (40) 

Moreover, when f is continuous, R + WR + is compact. 



Proof. The condition / = g gives immediately W (e l6 , e 1 ^) = W(e i< t' ,e l9 ), and so is self-adjoint. It can 
easily be seen that the Riesz projection R+ has distributional kernel 1/(1 — e l ( e ~^), and so W decomposes 
as 

W = M g [M f , R+] - M f [M g ,R + ] , (41) 
where all the operators are bounded. A simple calculation now shows that 

R+WR+ = (T gf - T g T f ) - (T fg - T f T g ), (42) 

and we apply the standard formulae Thk — ThTk = ^h(M)^k(z) an d = r^ 2 ) (see O p. 253]) to get equation 
(|40|) . The last statement follows by Hartman's theorem: the Hankel operators on the right-hand side of (j40| 
are compact when / is continuous. ■ 

Remark 3.2 We continue functions f £ L 2 to harmonic functions on ID) by means of the Poisson kernel, 
as in JSl p. 718]. 

Proposition 3.3 Suppose f = g G L°° , where g is holomorphic inside P. Then 

R+WR+=T* f T s . (43) 
Further, if R + WR+ has finite rank, then f is rational. 



Proof. Take / = g in Lemma T3. II to obtain the first part of the result. For the second part, note that 

Range(i?+l¥i?+) = Ker^T/^ = Ker(r / ) i = Range(r}) = Range(r /(2) ), (44) 

and apply Kronecker's theorem: Tk has finite rank if and only if k is rational, so ^ f(z) nas finite rank if and 
only if f(z) is rational, which implies that / is rational. ■ 
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